Entropy and Exact Matrix Product Representation of the Laughlin Wave Function 



O 

o 

(N 

Oh' 
(D 

cn 



S. Iblisdir\ J. I. Latorrei and R. Oms^ '^ 
^ Dept. d'Estructura i Constituents de la Materia, Universitat de Barcelona, 647 Diagonal, 08028 Barcelona, Spain 
^ School of Physical Sciences, University of Queensland, QLD 4072, Australia 
(Dated: 6th February 2008) 

An analytical expression for the von Neumann entropy of the Laughlin wave function is obtained for any 
possible bipartition between the particles described by this wave function, for filling fraction v = 1. Also, for 
filling fraction v = 1/m, where m is an odd integer, an upper bound on this entropy is exhibited. These results 
yield a bound on the smallest possible size of the matrices for an exact representation of the Laughlin ansatz 
in terms of a matrix product state. An analytical matrix product state representation of this state is proposed in 
terms of representations of the Clifford algebra. For u = 1, this representation is shown to be asymptotically 
optimal in the limit of a large number of particles. 

PACS numbers: 03.75.Ss, 03.75.Lm, 03.75.Kk 
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The puzzling behaviour of a cold gas of electrons confined 
in two spatial dimensions, and subjected to the effect of trans- 
verse magnetic fields, is the origin of two of the most stud- 
ied phenomena in condensed matter physics: the integer and 
the fractional quantum Hall effects While it is possi- 

ble to give a satisfactory explanation for the properties of the 
integer quantum Hall effect with filling fraction v = 1, by 
means of transport and related models pi, a complete under- 
standing of the fractional case is still missing. It is commonly 
believed that the interactions between the particles are essen- 
tially responsible for this situation. A relevant approach, in 
this respect, is the Laughlin ansatz for the wave function of the 
ground state of the system |:4]. While so far this state has only 
been proven to be an exact eigenstate of very specific Hamil- 
tonians @] and for some specific values of the filling fraction, 
it contains the relevant properties that the ground state of the 
real system must have. The LWF has also been turned to be 
a valuable ansatz in describing the physics of rapidly rotat- 
ing small atomic Bose-Einstein condensates 01. However, to 
check that the Laughlin wave function (LWF) is indeed the 
ground state of the fractional quantum Hall effect seems to 
be a difficult (computational) problem that has, so far, only 
been solved for systems made of a small number of particles, 
typically of the order of ten jsl. 

In this paper, we wish to provide some clues on why this 
problem is so difficult, in the light of recent developments 
in the fields of quantum information science and many-body 
physics. It is now known that a profound relation exists be- 
tween the difficulty to simulate numerically a quantum system 
in a given state and the von Neumann entropy between parts 
of the system this state exhibits, which can be understood in 
terms of the so-called Matrix Product States (MPS) 0]. With 
this relation in mind, we have computed the bipartite von Neu- 
mann entropy of a system as described by the Laughlin ansatz 
(i/ — 1). As we shall discuss, this computation is instructive 
in understanding the nature of the correlations exhibited by a 
system in such a state. To our knowledge, the behaviour of 
the von Neumann entropy was known only for the case of one 
particle versus the rest and in very specific situations 0]. In 
addition, we will provide an MPS representation for this state, 
which is optimal in the thermodynamical limit. This represen- 
tation will then be extended to some non-integer values of the 



filling fraction u. 

The LWF for 7i particles and filling fraction reads | 



. . . , 2„) = K(n) Y[{z, - Zj^/'Yle- 
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, (1) 



where Zj, j = 1 , . . . , n stands for the position of the particle 
j in the x — y plane written as a single complex coordinate 
Zj = Xj +iyj ■ Note that, but in a few cases, the computation of 
the normalisation constant, 'Ni, (n) seems to be a very difficult 
problem |9]. 

It is convenient to introduce an orthonormal monoparticle 
basis defined as 
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a = 0, 



, n — 1. 
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Using this basis, the LWF, for 1^=1, can be expressed as a 
Slater determinant: 



^'i(zi,...,z„) 
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°"0ai(zi) ■ • ■ 0a„(^;ri), 

(3) 

where e denotes the Levi-Civita completely antisymmetric 
tensor in n dimensions, with the convention £0>i^ - >"-i = 

The above expressions show that it is possible to correctly 
describe the quantum state by means of a local Hilbert space 
of dimension n for each particle. Since we have n particles, 
the total dimension of the complete Hilbert space scales then 
as ri". This superexponential scaling is at the heart of the 
difficulty to validate the ansatz Q by means of an exact diag- 
onalisation 

Let us now compute the von Neumann entropy for any sub- 
set of k particles for a system of n indistinguishable particles 
in the state . . . , 2„). Note that this von Neumann en- 

tropy cannot be interpreted as the number of distillable EPR 
pairs. Due to the symmetrisation, it is impossible to associate 
a label with the particles and perform the appropriate distilla- 
tion operations. 

The reduced density matrix for a subset of k particles out 
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of n reads 

Pk,n = Pk,n{wi, . . . ,Wk; Zi, . . . , Zk) = 
J dzk+1 ■ ■ ■ dZn . . . , Wfe, Zfe+1, . . . , Z„) X 

^l{zi, . . . ,Zk,Zk+l, ■ ■ ■ ,Zn). (4) 

From the orthonormality of the monoparticle basis, we get 

Pk,n — „! 2^a,6,c ^ X 

</'ai(w'l)'/'fcl(^l) • ■ ■K,iWk)(l)bA^k)- (5) 

This density matrix can be written in a diagonal form by in- 
troducing the following set of orthonormal basis for the subset 
of k particles: 

<^,iz) = -l^e-^---'-^--'-q>aAzi) . ■ ■ 4'a,{zk), (6) 

where (z) = (zi, . . . , z^) and where the indices are sorted 
such that Ck+i < ... < c„ and,thus, the combined index c 
ranges from 1 to ( ) . It is then possible to see that all eigen- 
values of pk n are identical, that is 



Pk^n - y7^^$:(w)$c(^). 

\k) c=l 

The von Neumann entropy then reads 

Sk^n = -Tr {pk^n l0g2 Pk^n) = loga 



(7) 



(8) 



and gets its maximum value for a bipartition of the system into 
two pieces of equal number of particles (this can be easily seen 
using strong subadditivity): 



Sk.n !i "S's. 



1 nn 



(9) 



Eq.(|9j shows that the effective dimension of the Hilbert space 
for half a Laughlin gas is Of2"), which is a weaker scaling 
than the naive growth, 0(n"'^), of the Hilbert space for n/2 
particles. 

Let us now perform a similiar computation in the case 
V ~ \ /to, when to = 2.s + 1 is an odd positive integer dif- 
ferent from 1. As explained in Reff^l. ^v{zi, ■ ■ ■ , Zn) can 
be expanded in terms of mutually orthogonal Slater determi- 
nants: 



Let A{n, s) denote the number of coefficients in this expan- 
sion. It is easy to see that the rank Tfc of the reduced density 
matrix of any subset of k particles is bounded as 



Tfc < A{n, s) 



(11) 



Indeed, as we have seen in the study of the case v = \, (JJ) 
is the Schmidt rank of each of the individual Slater determi- 
nants, and the Schmidt rank of a sum of kets is smaller than 
or equal to the sum of the ranks of individual kets. As abun- 
dantly discussed in Ref.lIP], to give a closed form for A{n, s) 
seems to be a very hard problem. Fortunately, an upper bound 
on A{n, s) useful for our purposes can be derived as follows. 
Due to the constraints between the indices, the expansion ([TOJ 
doesn't feature more than ((2s + l)(n — 1) + 1)" terms. This 
trivial bound could be reached, would each index be allowed 
to range from to (2s + l)(n— 1) independently from the val- 
ues taken by the other indices. But one can further constraint 
A(n, s): since a given permuation of the columns of a deter- 
minant produces the same determinant (up to a sign factor), 
we have: 



A{n, s) < ((2s + l)(n - 1) + l)"/n! 
Therefore, for k = n/2 and in the large n limit, 

(2s + 



(12) 



Tn/2 < 



2" = (4s + 2)", (13) 



and the von Neumann entropy obeys 

5^,„ <0(nlog(2.s + l)). 



(14) 



This bound can actually be slightly improved upon using the 
fact that all indices Iq, . . . , Z„_i assume different values 

The dimensions of the matrices in an MPS representation 
of a given quantum state are related to the von Neumann en- 
tropy of bipartitions of the system when its degrees of freedom 
are ordered on a line 0]. According to our previous calcula- 
tion, an MPS representation of the LWF in terms of matrices 
of size 0(2") should be possible. More precisely, it should 
be possible to find a representation of the coefficients of the 
LWF written in the monoparticle basis in terms of products of 
matrices: 

n-l 

*,(zi...z„) = — J2 tr(AW''i ...AW'^") X 

ai , . ..,an—0 



^^{zi...z„) = Y[e 



(^) 



h,...,i„ 



(10) 



where the indices /i, . . . , /„ are constrained by < /i < 
. . . < Z„ < (2s + l){n - 1). 



(15) 



where the matrix is associated with the particle i be- 

ing in the monoparticle state (^ai ■ These matrices have a size 
which is the same for all values of i, a^: x ^ X- We will now 
see that the properties of Clifford algebras are well suited in 
order to find these matrices. The Clifford algebra CZ(0, n) is 
defined by 



{7^7''} = 25 



ah 



0, 



1, 



(16) 
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where each matrix 7°^ has indices a, f3 = 1, . . . ,x (see for 
example fl^ l. 

Let us start with the case where n is even. The representa- 
tion theory of the CHfford algebra dictates that x ~ 2"/'^. The 
matrices 7°^ provide the following MPS construction: 



(2i)"/- 



,tr(7''^...7''"75). 



(17) 



j)"/2 -yO . . . --y" 1. This result emerges from 



where 75 e 

the basic trace properties of the matrices 75 and 7°, a = 
0, . . . ,71 — 1, and shows that all the matrices A^'l"' can be 
taken the same for all particles but one. For example: 



EE 7"^' i = 1, , 



,71—1, 



(18) 



Our construction shows that, in spite of the symmetry of the 
system, the MPS representation is not made of a set of identi- 
cal matrices. This property was also observed for W-states in 

m. 

The MPS construction ( I18t is asymptotically optimal in the 
sense that the matrices have, in the limit of large n, the mini- 
mal size. To begin with, let us first compute how large should 
the matrices be in an open boundary MPS representation. The 
rank of the indices of the different matrices then corresponds 
to Schmidt numbers of bipartite decompositions of the quan- 
tum state 1 7]. Given a bipartition of the system into two con- 
tiguous pieces, and assuming that the x eigenvalues of the 
corresponding reduced density matrices of the two subsys- 
tems are all equal to 1 /x, we get a relation between x and the 
maximum possible von Neumann entropy for that bipartition: 
Sm-dK ~ log2 X- In general, though, the eigenvalues may not be 
equal, and therefore x ^ '^^ (open boundary MPS), S 
denoting the actual von Neumann entropy of the bipartition. 
This argument fails for periodic boundary MPS since corre- 
lations between two sets of a bipartition can flow across two 
boundaries. Yet, it is possible to take the periodic-boundary 
MPS as a product of two effective matrices L^pRpa (s.g. for 
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a bipartition of the system Lq,/ 
and an analogous construction for the right part of the sys- 
tem represented by i?^ct). A simple Schmidt decomposition 
of the LapRpa construction shows that the Schmidt rank will 
now run up to y^. The same argument as for open boundary 
conditions now shows that 



X 



periodic boundary MPS. 



(19) 



The above considerations confirm that, in the case of pe- 
riodic boundary conditions, it is equivalent to argue that the 
effective rank of a bipartition has contributions coming from 
the two borders. As we have already seen, in our case the 
maximum possible entropy over all bipartitions corresponds 
to S !i ^ n, which in turn implies that the dimension of the 
matrices of the Clifford algebra precisely matches the lower 
bound provided by this entropy for large n 



X 



2ri/2 



(20) 



Therefore, the representation (I18t is optimal in the limit of 
large n values. 

Explicit representation of the MPS construction can be ob- 
tained from the chain of isomorphisms 



CZ(0, 71 + 2) « CZ(0, 2) ® Cl{Q, 71). 



(21) 



For 71 = 2, set 7*^ = a^, 7^ = and 75 = a^. A represen- 
tation for 71 + 2 is constructed from a representation of 71 as 
follows: 

7(„+2) = 1®7(„), i = 0,...,n-l, 

7(„+2) = 7("2)"'^75,(«), i = 71, 71 + 1. (22) 

A representation of the Clifford algebra for the case 71 odd 
can be simply derived from a representation for 71 — 1, upon 
taking 



7(„) = 7(„-i), i = 0, ...,7i-2, 

n— 1 

7(„) ~ 75(n-l)- 



(23) 



Then, it can be shown that 75(„) — (— and that the 
antisymmetric tensor is given by 

_jl/2 

ai...a„ trf-T^i ^°") (24) 

(2i)("-i)/2 ^ ' 

It is also possible to extend the MPS construction of the 
Laughlin wave function we have obtained for the case v ~ 1 
to the fractional case whenever v ~ where m is an integer. 
The basic idea is to use the property of the product of traces 
Tr(Ai . . . A")Tr(Bi . . . B") = Ty{{A^®B'^) . . . (A"(g)B")) 
in a recursive way. Let us exemplify the construction in the 
case of 1/zy = 2, (bosonic statistics). 

Up to a global prefactor and a set of particle-dependent fac- 
tors, the problem of finding a faithful MPS representation es- 
sentially reduces to representing 



tr 



(7''"...7"^"75)tr(7"5r..7'^^N5)x 



The structure of the coefficients is then correctly repre- 
sented by 



where, for all z. 



(26) 



®2 



ai+c[2 — c 

(27) 

It follows that the dimensions of the combined matrices pro- 
duces an effective x = 2". In general, an MPS representation 
of the LWF with filling fraction i/ = 1 /m for integer 7?i can 

be defined in terms of the matrices P^/,„ = (7^) ™, 



^l/rn - E 



I^i±:::±^7'^^®---®7'^-, (28) 
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yielding an effective x = 2t^. This value of x is far above 
what one could expect to be the minimum from Eq. (ll4> so 
that a much more economical construction should exist. In- 
deed, in the case n = 2, we have found the following con- 
struction. The matrices F? , can be viewed as acting on 

1 / m 

the TO-fold tensor product representation of SU(2). Define 
'^i/m ~ where G is a projector onto the spin- 

m/2 irreducible subspace contained in (1/2)®'" (a piece of 
Clebsch-Gordan matrix). After some algebra, we find that 
the explicit form of each matrix has just a single entry 1 in 
its antidiagonal, T^jy^^(c + l,n — c) = 1, and the rest is 0. 

Further let T:^^^^ denote an rn + 1 x m + 1 matrix whose di- 
agonal elements are defined as {T^l^^)ii ~ ^iid 
whose off-diagonal elements are zero. We have numerically 
checked and yield an exact matrix product repre- 

sentation of in the case n ~ 2 for m ~ 1 ... 15, and we 
believe that they also do for arbitrary values of m. Note that 
for n — 2, the von Neumann entropy of the state of one par- 
ticle can be easily calculated for any value of = 1/m. One 
gets 5*1 /m = X] jLo (7) ^^^2 ( j*) ■ Clearly the size of the new 
matrices is now m + 1 instead of 2™, that is exponentially 
more economical. We should note that this representation has 
an interest as a starting point to get intuition on more involved 
cases. But, as such, it is not very useful in practice since the 
computation of a mean value (Vi (8) V2) requires an effort that 
scales as O(m^), which is the same as if using the Schmidt 
decomposition for '9^. 

For some filling fractions ly, the foregoing analysis can 
be easily extended to study excited states of systems whose 
ground state is described by an LWF Let us consider the sim- 
plest case of one quasiparticle localised at a position za- The 
corresponding wave function reads ||6| 

n 

*[f^l(zi,...,z„) =[](z,-z^)«',(zi,...,z„), (29) 

i=l 



Now since 

j_2 1 ■ ■ ■ ^" 

we have that ^^^I'^^zi, . . . , Zn) = 
e'^^''-'"z]Czl^ ...zl^"^^{zi,...,Zn),v/ithiy' = vjiy+l). 
With calculations akin to the ones performed above, an 
MPS representation of is easily derived. Defining 

Y'^[za) = jy^^=o l^^^lt ^ have (up to normalisation): 

*|;,^i(zi,...,^„) = tr(r0(z^)r'^^..r°"r5) 

(j>ai{zi) . . .(I)a„iz„). (30) 

MPS representations of m-quasiparticle excited states can be 
computed likewise. 

In summary, we have computed the von Neumann entropy 
of a Laughlin gas for any bipartition {k,n ~ k) of the system 
under study and for various values of the filling fraction v. We 
have seen that this entropy grows at most linearly with n and 
logarithmically with i^. Since this quantity may be related to 
the difficulty of numerically simulate the gas, this computa- 
tion sheds new light on why such systems are so difficult to 
study. Next, we have provided an MPS representation of the 
Laughlin wave function, that is asymptotically optimal in the 
case — 1. We believe this representation can be exploited 
to compute various quantities related to the LWF, such as its 
norm. This will be the subject of further investigation. 
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